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We study cubic ferroelectrics films that become uniaxial with a polar axis perpendicular to the 
film because of a misfit strain due to a substrate. The main present result is the analytical account 
«j ' for the elastic anisotropy as well as the anisotropy of the electrostriction. They define, in particular, 

, an orientation of the domain boundaries and stabilizing or destabilizing effect of inhomogeneous 

elastic strains on the single domain state. We apply the general results to perovskite systems like 
BaTiOa/SrRuOs/SrTiOs films and find that at least not far from the ferroelectric phase transition 
the equilibrium domain structure consists of the stripes along the cubic axes or at 45 degrees to 
them. We have also showed that in this system the inhomogeneous strains increase stability with 
regards to the small fluctuations of the metastable single domain state, which may exist not very 
close to the ferroelectric transition. The latter analytical result is in qualitative agreement with 
the numerical result by Pertsev and Kohlstedt [Phys. Rev. Lett. 98, 257603 (2007)], but we show 
that the effect is much smaller than those authors claim. We have found also that under certain 
conditions on the material constants, which are not satisfied in the perovskites but are not forbidden 
I either, a checkerboard domain structure can be realized instead of the stripe-like one and that the 

^ ■ polarization-strain coupling decreases stability of a single domain state instead of increasing it. 

The single domain state is metastable at certain large thicknesses and becomes suitable for memory 
applications at even larger thicknesses when the lifetime of the metastable state becomes sufficiently 
large. 
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PACS numbers: 77.80. Dj, 77.80.bn, 77.55.Px 



Ch ' I. INTRODUCTION 

o ■ 
o 

Properties of domain structures in thin ferroelectric films is currently a focus of extensive research. It is expected, 
quite naturally, that an understanding and an ability to control these properties will determine the prospects of 
applications of nanometer-size ferroelectrics. It depends critically on the external conditions like presence or absence 
■ of electrodes. In this paper, we discuss domain structures in a system, which is, perhaps, the most important for 
iy~j I applications: a ferroelectric film with electrodes. The polar axis of the material is perpendicular to the film plane 
. and the electrodes are 'real' meaning that the electric field penetrates into them, although only over tiny depths 
t-H ' < 1A. This is an adequate model for the perovskite ferroelectric films on a substrate with compressive strain, like 
! BaTi0 3 /SrRu0 3 /SrTi0 3 (BTO/SRO/STO)^ where the misfit strain drives the FE film into a uniaxial state. We 
supplement our analytical results with the relevant numerical estimates for BaTi0 3 (BTO), PbTi0 3 (PTO), and 
t-H , Pb(Zi'o.5Tio.5)0 3 (PZT) using the material constants available in the literature 

Incomplete screening of the depolarizing field by SrRu0 3 electrode leads to an absolute instability of a single domain 
state and formation of a sinusoidal domain structure when thickness of BaTi0 3 film is slightly above the minimal 
thickness compatible with ferroelectricity in this system^—. It seems that this situation is typical of real electrodes 
and we shall consider this case only. To find the minimal thickness, one does not need to take into account higher 
order terms in the Landau- Ginzburg-Devonshire (LGD) free energy including the terms describing the electrostriction 
since the problem of stability of the paraelectric phase is linear—. But in order to reveal the characteristics of this 
structure, i.e. to find out if the equilibrium structure is stripe-like or checkerboard and how the domain boundaries 
are oriented one has to take into account the anizotropy of elastic and electrostrictive properties of the ferroelectric. 
This is the main goal of the present paper. Specifically, we consider the case of cubic crystal anisotropy of clastic 
and electrostrictive properties only. This is relevant for films of cubic perovskites which become tetragonal because of 
in-plane misfit strains due to cubic substrates like in the above-mentioned system. The change of cubic anisotropy to 
tetragonal affects most strongly the dielectric properties since the crystals are "soft" dielectrically. They have much 
smaller effect on the elastic and electrostrictive properties, which can be considered to be the same as in cubic parent 
crystals. 

Explicit account for the electrostriction and anisotropic elasticity is relevant also for study of stability of single 
domain state. This has been correctly pointed out by Pertsev and Kohlstedt£«i although these authors have missed 
several important points. Importantly, however—, they overlooked that the state whose stability they were studying 



2 



was actually metastable. Therefore, its stability with respect to small fluctuations did not mean that this state 
can be used in memory applications. Indeed, its lifetime is very short if the film thickness is not sufficiently larger 
than that calculated by Pertsev and Kohlstedt as the limit of single domain stability. We shall also discuss this 
stability among other questions This makes sense because of several reasons. First, Pertsev and Kohlstedt performed 
numerical calculations using material constants for BaTiC>3 and Pb(Zro. 5^0.5)03 and the electrode parameters of 
SrRu03 while our results are analytical and apply to any material of the same symmetry. Moreover, our method 
applies to other symmetries as well. Second, Pertsev and Kohlstedt studied stability of the single domain state with 
respect to "polarization wave" -like fluctuations with a single specific direction of the wave vector, while we consider 
waves with fc-vectors in arbitrary direction. Third, Pertsev and Kohlstedt apparently misinterpreted their own results 
by mixing together the well-known effect of homogeneous misfit strains and the effect of strains due to inhomogeneous 
polarization. In fact, the misfit strain simply results in renormalization of the materials constants and was effectively 
taken into account by all the previous authors. Only the account of the inhomogeneous polarization was pioneered in 
Ref,i2. We show that this effect was vastly overestimated by Pertsev and Kohlstedt. In fact, the formal difference by 
more than an order of magnitude between the results with and without account for electrostriction that they claimed 
stems from improper comparison of the compressed film with the materials constants renormalized by the misfit strain 
to one without any such renormalization at all, and not from the effect of the inhomogeneous strains on stability of 
single domain state. 

Our analytical calculations provide a general view on the role of the inhomogeneous strains in stability of single 
domain state. In particular, they reveal a possibility which seems academic at the moment but no reason is seen 
to exclude it altogether. We mean a specific state where elasticity provokes domain formation of non-ferroelastic 
ferroelectric 180° domains. Such a state is realized if a certain condition on the electrostrictive and elastic constants 
is met. We are not aware of an experimental realization of these conditions but we cannot find arguments prohibiting 
them. It is worth mentioning that a qualitative conclusion about possibility of both stabilizing and destabilizing 
role of inhomogeneous elastic strains for single domain state in ferroelectric films on substrates has been made in 
our previous paper where we considered an academic case of a single electrostriction constant and assumed isotropic 
elasticity^. A surprising result of the present work is that the destabilizing effect of the inhomogeneous strains may 
be very large contrary to the stabilizing one. Another unexpected result is the possibility of a checkerboard domain 
state if some conditions on the material constants are met. Let us mention that without account for the anisotropic 
polarization-strain coupling one comes to the conclusion about impossibility of such a state. For the perovskites this 
conclusion remains valid but not in the general case. 

Studying the sinusoidal domain structure in BTO, PTO and PZT films on SrTi03,we find that the equilibrium 
orientation of the "domain walls" is parallel (perpendicular) to the cubic axes in the film plane for BTO and PTO and 
is at 45° to these axes for PZT. In all cases the free energy of the sinusoidal domain structure depends very weakly on 
the domain wall orientation. This is mainly due to both systems being nearly isotropic elastically and, additionally, the 
relevant electrostriction constant is relatively small. This observation may be important for understanding domain 
creation at smallest thicknesses of the ferroelectric films. For BTO/SRO/STO system, our analytical calculations 
provide confirmation of the qualitative result of Pertsev and Kohlstedt about stabilizing effect of inhomogeneous 
elastic strains for single domain state in this system but with the above mentioned strong disagreement with their 
statement about importance of this effect. 

Having mentioned advantages and new possibilities provided by analytical calculations we should mention also their 
inherent shortcomings. Our analytical method is feasible within a certain approximation only. This approximation 
implies that the domain period is less than the film thickness. This condition is fulfilled for thick enough films but 
in very thin films the two quantities are in fact comparable. Therefore, the accuracy of our calculations should be 
investigated for these films, so the new numerical studies are desirable. We do not expect, however, that the difference 
between the results of approximated and more exact calculations either within a continuous medium theory or within 
microscopic theories will be very large given close results of continuous and first principles theories even for films that 
are just several unit cell thick (see, e.g.^). 

The paper is organized as follows: we describe the approximations used and define the terms in the LGD free energy 
that can be neglected within our approximation in Sec|TTJ This let us avoid unnecessary lengthy formulas in the rest of 
the paper. We spell out the constituent equations in Scc lIIII and then solve the general problem for the 'polarization 
wave' (embryonic stripe domains) in the FE film with full account for elastic coupling. This is further used in SeclVl 
to determine that the domain walls align with crystallographic cubic axes. Then, we find the conditions when the 
monodomain state first loses its stability with regards to the stripe domain structure in Sec lVIl One previously 
unexplored possibility is that the system can lose stability with regards to checkerboard domain structure, but our 
results in Sec lVIII show that such a structure is absolutely unstablein perovskites although it is not necessarily so in 
the general case. We summarize the present results in the Conclusions. 
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FIG. 1: Schematic of the (perovskite) ferroelectric film with thickness I and metal electrodes (with screening length A) on a 
misfit substrate. The misfit makes the film a uniaxial ferroelectric with a spontaneous polarization along z— axis. 



II. OUTLINE OF THE METHOD AND THE APPROXIMATIONS USED 



The main conclusions of this paper are made by analyzing the formula for free energy of the total system as a 
function of the amplitude a of the ferroelectric "polarization wave" presenting the sinusoidal domain structure and 
the homogeneous part of the ferroelectric polarization, p. For the electrode and the film parameters of a system like 
BTO/SRO/STO, the ferroelectric polarization that is perpendicular to the film plane, schematic of which is shown 
in FigQ] has the form 



P z (x, y, z) = p + a cos kr cos qz, 



(1) 



where the orientation of k in the x, y plane is not fixed, 27r/fc is the period of the sinusoidal domain structure, q = it /I, 
and I is the film thickness. To find the desired free energy, F(a,p), one has to find the elastic strains and the non- 
ferroelectric polarization Pi = (P x , P y ) as functions of a and p to present the total free energy as a function of a and 
p only. The total free energy contains contributions of the ferroelectric film, of the substrate and of the electrode. 
In principle, it should also contain a contribution of the voltage source but we consider here a short-circuited system 
and are not concerned with the latter contribution. 

When calculating elastic strains in the ferroelectric which accompany the inhomogeneous polarization forming the 
sinusoidal domain structure, we follow the same philosophy as in our previous workJ^. In principle, when calculating 
these strains the inhomogeneous strains in the substrate should be taken into account. However, it is well known 
that they propagate into the substrate for about the same distances as the scale of inhomogeneity in the film (x, y) 
plane. In our case, these inhomogeneities are due to the domain structure, i.e. this scale is the period of the domain 
structure. Then, it is convenient to consider relatively thick films since the period of the domain structure is relatively 
small, specifically, it is much less than the film thickness^, Fig. [2] The contribution of the substrate is its elastic 
energy which, as we have mentioned above, is concentrated within a volume which is much smaller than the film 
volume, as defined by a small factor q/k = ir/kl <C 1. Another convenience of the thick film limit is that it is possible 
to disregard the boundary conditions for the inhomogeneous parts of elastic strains and stresses at the surfaces of the 
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FIG. 2: (color online) Schematic of the ferroelectric film on the misfit substrate at the onset of sinusoidal polarization wave. 
The elastic coupling to the substrate allows inhomogeneous deformations but prohibits homogeneous strains in plane of the 
film. 



ferroelectric. Indeed, if we obtain a solution, which does not satisfy the boundary conditions, we can find corrections 
to such a solution in a way that is conventionally used in the elasticity theory. First, we apply the external forces 
to the surfaces, which are necessary to meet the boundary conditions with the strains corresponding to our solution 
making this solution correct. Second, we apply forces opposite to the previous ones and find the strains produced 
by the new forces. These strains provide the correction to the original solution we were looking for. Once more, it 
is sufficient to observe that in our case the external forces have the period of the domain structure to understand 
that the elastic energy associated with the corrections necessary to satisfy the boundary conditions can be neglected 
quite similarly to the elastic energy of the substrate. Another convenience of thick film approximation is given by the 
possibility to neglect those terms in the LGD free energy, which describe the electrostriction but contain components 
of non-ferroelectric polarization. According to Refs.— &, 



_P_l (x, y, z) = (k/k)a± sinfer sinqz, 



(2) 



where a± w aq/k, q/k = n/kl <C l.The electrostriction terms in the LGD free energy with non-ferroelectric compo- 
nents of polarization may contain the ferroelectric component, like P x P z u xz , or may not contain them, like in the 
term P x P y u xy . In both cases, they contribute to a 4 and p 2 a 2 terms in the free energy depending on a and p. In the 



first case, this contribution is proportional to (q/k) and in the second to (q/k) . Since there are also the terms 
p 2 a 2 that do not contain the small factor q/k, the contribution of these terms can be neglected. 
Taking this into account, we write down the LGD free energy in the form: 



F (P, u lk ) = Fx (P) + F 2 (u lk ) + F 3 (P, u. 



ik ) 



where 



Fl (p) = ^p 2 + ^+ l -G(v^p z y 



-kP, 



A 



bz 



± P 2 



(3) 



(4) 



F 2 (u ik ) = -Ai (u 2 xx + u 2 yy + u 2 zz ) + A 2 (u xx u yy + u xx u zz + u zz u vv ) + 1[i (u 2 xy + u 2 zy + u 2 xz ) 



(5) 



F 3 (P, u lk ) = qnu zz P 2 + q 12 (u xx 



l) ± z 



(6) 



Here, (711(12) are the standard piezo-electric coefficients that should not be confused with the parameter q defining the 

transversal profile of the polarization wave P. In Eq.g]), A = j(T-T c ),B,G = const, ^ j_= (d x , d y ) the gradient 
in the plane of the film, Pb z is the non-ferroelectric ('base') part of the polarization perpendicular to the electrodes^, 
A± > 0. Following Refs^^, we have neglected a term with the gradient in z— direction since it is much smaller than 
the one in plane of the film, d z <C vj_. It is worth mentioning that we have not included the energy of the electric 
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field. into the LGD free energy. The reason is that we shall use it to write down the constituent equations only. 
We shall eliminate u^. , P±, P bz as well the electric field components from the system of constituent electrostatics 
equations to obtain two coupled equations of state for a and p. We shall obtain F (p, a) from the resulting equations. 
This is possible because of the thick films approximation. The most straightforward method to obtain F (p, a) would 
be to substitute Eq.([T]) into Eq.Q supplemented by the electric field energy and to integrate over the film volume. In 
general, the result would not be the same as the one obtained from the constituent equations because of approximate 
character of Eq.([T]). However, for q — ir/l the two results coincide and that makes it possible to use a more convenient 
method of the constituent equations. 

III. CONSTITUENT EQUATIONS 

For the polarization components one has: 

AP Z + BP* - GV\P, + 2q n P z u zz + 2q 12 P z (u xx + u yy ) = E z , (7) 

Pbz = kE z , (8) 
P x = A X E X . (9) 

Before writing down the equations for the strain, we shall eliminate the electric field from the above three equations. 
Assuming Eq.([T|) for P z , Eq.@ for P^ and putting^: 

Eqz = E + E z cosfcr cosqz, E± = (k/k)E k ± sin kr sin qz (10) 

we can replace Eqs.©,© with 

Ap + [BP* + 2q n P z u zz + 2q 12 P z u ±± ] hom = E , (11) 

(A + Gk 2 ) a + [BP* + 2 qil P z u zz + 2q 12 P z uxi_] cc = E h z , (12) 

A ±a± =Ei, (13) 

where u±± = u xx +u yv , [. . -] hom and [. . .] cc denote the homogeneous part (k = 0) and the part proportional to 
cos kr cos qz of the expression in the brackets, correspondingly. Of course, as a result of this replacement, a part of 
the l.h.s. of Eq.Q is lost but it corresponds to the higher harmonics of the sinusoidal distribution of the polarization 
and these harmonics can be neglected close to the transition^. 

The homogeneous part of the electric field E z0 can be calculated as, e.g., in ReL 6 yielding for the short-circuited 
case 

4ird , , 

E 0z = j— -:p, 14 

where d is the thickness of the dead layer and e e its dielectric constant. Recall that real electrodes have finite albeit 
small (Thomas-Fermi) screening length A, which is completely analogous^ to a presence of the 'dead' non-ferroelectric 
layers at the interface with thickness d/2 = A. Using Eq. (TUT) , the equation (fTTj) gets the form: 

A x p + [BP* + 2 qil P z u zz + 2q 12 P z (u xx + u yy )] hom = 0. (15) 

where 

4ird . And „. 

Ai = A+ ^A + —, (16) 

e b d + e e l e e l 



since usually the dead layer is very thin, Ebd <C e e l. To transform Eqs. (II 2|) we use the electrostatics equation, 

divD = 0, (17) 

where D is the dielectric displacement, firstly for the ferroelectric material, taking into account that D = 
(e±E±, £bE z + 4irP z ), where e± = 1 + 4:Tr/A±, and — 1 + 4tt/k is the 'base' non-critical dielectric constant&il, 
and together with the equation curl_E = 0, we find that 
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e±k^ 

Substituting (jl8p into the equation for the amplitude of the 'polarization wave' a, Eq. ([T2|) . we rewrite the latter as: 

A + Gk 2 + a + [BP* + 2q n P z u zz + 2q 12 P z u ±± ] cc = 0. (19) 

The nontrivial solution of the above equation first appears when the coefficient in the first term in round brackets 
before the amplitude a in the above equation first crosses zero, i.e. when A = [-Gk 2 - Airq 2 / (e±k 2 )] max . That 
takes place at some A < 0, so upon lowering temperature at constant thickness the transition into domain state 
occurs somewhat below the bulk critical temperature T c , in other words. The transition for varying thickness of the 
film at some constant temperature T < T c takes place when the thickness exceeds some critical value. Therefore, 
the first nontrivial solution appears for the 'polarization wave'^- with the wave number k that minimizes the sum 
Gk 2 + Airq 2 / (e±fc 2 ), so that (recall that q = n/l) 

-(Sf "(i&f ■ 

We can now rewrite Eq. (|12p as the homogeneous one: 

A 2 a + [BP? + 2q n P z u zz + 2q 12 P z u ±± ] cc = 0, (21) 

where 

A 2 =A + 2Gk 2 . (22) 
It is seen from Eq. (J6j) that the only source of elastic stresses and strains is P 2 (x, y, z) in our approximation. Since 

a 2 

P 2 = p 2 + 2pa cos kr cos qz + — (1 + cos 2qz + cos 2kr + cos 2kr cos 2qz) , (23) 
we should expect that 

u zz — uf} + u^J cos 2qz + u^ 2 ) cos kr cos qz + u^J cos 2kr + uyj cos 2kr cos 2qz, (24) 

while for u XXl u yy , and for u±±_ we shall have similar formulas with the homogeneous part (first term in the above 
expression) absent because of the substrate. The superscripts (0) — (4) denote contributions with different types of 
the coordinate dependencies as defined by Eq. (j2~4"|) . Below, we use the same superscripts for both the coefficients and 
the functions. 

Substituting Eqs. (f2"3|) , (j24l and analogous equations for u xx and u yy into Eqs.([T5 l) . ([2T]) . we find: 

A lP +[BP?] hom + 2 qil (puW + a -f-\ +« 12 |^i 2 i = 0, (25) 



A 2ffl +[BP z 3 ] cc + 2 9ll 



pu<- 2 ] + a I «w + 



(o) 



(7) 

We shall calculate the values u/ k in the next Section by solving the elastic problem explicitly. 

Importantly, the above equation of state ([25]) suggests that the film would tend to transform into a single domain 
(SD) state with p 7^ and a = at temperature T^ D such that A\ — or, in other words, 



(1) , (3) 



(4)~ 



■2qi 



' , , /7, (1) +7, (3) 7, (4) ' 



= 0. (26) 



A (Tg ) = —And/ (e e l) . 



(27) 



The second equation of state (|26|) yields a transition into a domain state (p = and a ^ 0) at the temperature Td 
such that A 2 = 0. or 



A{T d ) = -2Gk 2 = -4 



5 G 



1/2 



1 



(28) 
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Recall that in the present case, corresponding to BaTiOs/SrRuOa/SrTiOa 5 , 

l -">2Gk*~^, (29) 



■i eTi 



where d at = \ZnG ~ 1A is the small 'atomic' length scale (G = 0.3A 2 for BaTiC>3 5 ' 6 ). The above relation means that 
the paraphase gives way to the domain phase, with a / 0, thus preventing it from reaching the temperature where 
it could have transformed into a single domain state. Obviously, same is true of the phase transformations in the film 
as a function of thickness at constant temperature. There, one can introduce the critical thickness for domains, Id, 
where 

7T 3 G\ 1/2 1 



A = -2Gk z = -4 -, (30) 

V e -L / Id 

and the 'critical thickness for the single domain state' lf D , such that 

A = -And/ (eJ s c D ) . (31) 
These introduced critical thicknesses and temperatures are discussed in detail below in Sec lVIl 



IV. ELASTIC PROBLEM 

Using Eqs. ([5]), ([5]), we obtain for the diagonal components of the elastic stress tensor: 

a xx = \iu xx + A 2 (u yy + u zz ) + q 12 P z , (32) 

<7yy = \\Uyy + A 2 {U XX + U zz ) + q^P^, (33) 

a zz = Xiu zz + A 2 {u xx + Uyy) + qi\P z , (34) 



and formulas of the type 



°~xy — '2^U X y 1 (35) 



for the off-diagonal components. 

We have already mentioned that the only itifc component which has a non-zero homogeneous part is u zz . This part 
is easily found from the condition at the free surface: <j zz = at z = 1/2. From Eq. (f3~4"]) . one finds: 

For the parts depending on z only, the equations of elastic equilibrium take the form: 

dvg>/dz = 0, (37) 

i.e. — const = since it should vanish at the free surface (z = 1/2). Therefore, Eqs. (|3~4")) , (|2"31 yield 

u« =- 9ll a 2 /(4A 1 ), (38) 

and 

= <J = °. ( 39 ) 

because of the Saint- Venant's elastic compatibility conditions for z— only dependent strains. 

When solving the rest of the elastic problem, we shall use the small parameter q/k <C 1. This allows us to neglect 
the derivatives with respect to z: formally, d/dz <§; d/dx, d/dy. As a result, the equations of the elastic equilibrium 
acquire the form: 

5(J (2-4) Sff (2-4) 9fT (2-4) 9(T (2-4) ^t 2 " 4 ) gj 2 ^ 
u "xx Uo X y _ UOyz OCTzx _ <J<Jyy OOyx _ ^ 

dx dy dy dx dy dx 
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where the superscripts (2 — 4) denote the part of stresses that is due to the three last terms in Eg. (l23|) . which we 
denote as P 2 ^ 2 ~ 4 ^. Explicitly, 

, d 2 4 2 ~ 4) . a 2 4 2 " 4) s 2 4 2 " 4) 9 2 4 2 ~ 4) 5p 2(2 - 4) n 

Ai 5-5 h A 2 „ „ h A* KM n o ^912 5 = °) ( 41 ) 

d 2 4 2 ~ 4) 9 2 4 2 - 4) n 
"^^ + '^^= < (42) 

. 9 2 4 2 ' 4) , » d 2 4 2 ~ 4) , dur i] , 9 2 4 2 ' 4) ^ 2(2 ' 4) n 
oy z oyox oxoy ox* ay 

Analogously to the isotropic case 1 ^, we shall put the conditions 4 2 ^ — that satisfy Eq.([J3]) but not, of course, the 
boundary conditions. The latter is not important in our approximation, as we argued above. Therefore, we conclude 
that 

ug> = «£> = i&> = 0, (44) 

and we are left with only two equations to solve. It is convenient to solve them separately for (2) and (3, 4) parts, 
since they correspond to different spatial harmonics. 

Simplifying the remaining equations ([4"T j) . (T4"3")) , we obtain: 

. a 2 ^ 2 ' 4) » ^ 2 4 2 ~ 4) ^ 2 4 2 - 4) ap 2(2 - 4) n 

Al ^~ + (A2 + ^-dyd^- + + qi2 ^x~ = °' (45) 

x 5 2 4 2 ~ 4) M ,d 2 4 2 ~ 4) a 2 4 2 ~ 4) ap 2(2 ~ 4) n 

Ai " +(A 2 + m) * +M " +gi2^ =0. (46) 

ojr oyox ox z oy 

For the terms u^ 2 \ we have p 2 ' 2 -* cx cosfer, d x ( y )P 2 cx —k x ( y j sinfcr, meaning that u x t y ) cx sinfcr, d 2 Ui/dx 2 — —k 2 Ui, 
etc. Then, 

Aifc 2 4 2) + (A 2 + fi)k x k v u y 2) + iik 2 y u x 2) + q 12 k x 2pa = 0, (47) 

Aifc 2 4 2) + (A 2 + H) k x k y u x 2) + iik 2 x u {2) + qi 2 k y 2pa = 0. (48) 

The terms u^'^ correspond to higher spatial harmonics in Eq. (f2U)l . but they should be taken into account since they 
contribute to the terms with the main harmonic in the constituent equation (|26[) . Since for this part P 2 cx cos2/cr, 
we obtain a slightly different set of equations: 

\ik 2 x u x 3 ^ + (A 2 + p)k x k y u y 3 ^ + M 2 4 M) + ?i2fcxy - 0, (49) 
A!fc 2 4 3 ' 4 ) + (A a + n) k x k y u x z ^ + /ifc 2 4 3 ' 4) + 9i 2 fcyy = 0. (50) 

Note that for the constituent equations we need the combinations: 

,.(2) _ ,.(2) , (2) _ , (2) , . (2) 

<4 3 i 4) = 4 3 / ) +4 3 / ) = 2 ^4 M) +2fc,4 M) - 

We find: 

ui 2 { = -2< Zl2 ap/ (0) , (51) 
where 9 is defined by k x = fccos#, k y = ksinO and the function / {&) is 
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f(0)=2- 



(Ai - A 2 ) sin 2 29 + 2/z cos 2 29 



' (Ai + A 2 + 2fi) (Ai - A 2 ) sin 2 29 + 4A x /x cos 2 26> ' 
For terms, corresponding to cos 2kr and cos 2qz, we obtain 



(52) 



„(3) _ „(4) _ a 

u ±± - u ±± - -qi2 — J (9) . 



Using the results of solution of the elastic problem, Eqs. (|3l^ . (j4"4")) . (|51[) . (fi33")) . we can write Eqs. ([23|) . (l2^1) as 



A lP + B 



Hi 



p A + pa 



B 



4 ( q \ 



3 V 2A 



J-+9?a/(f)) 



(53) 



(54) 



9 / 4 

A 2 a + a 3 — B 

16 V 3 



5ii , 5i2 
Ai 2 



/(*) 



ap z 3 B 



2Ai 



0. 



(55) 



From these two constituent equations corresponding to an extremum of the free energy, one can easily reconstruct 
the free energy: 



,-. , 4 5 A 2 2 B a 3B 1 (9) 2 , 9B 2 (6») 4 



where 



Ai o 



Ai 



1 2 i2f(0) 



, B 2 (9) = B 



r„2 
Ai 



(56) 



(57) 



are the Landau coefficients before the quartic terms renormalized by the strain. The form of the free energy is the 
same as in the isotropic case£&i£, but, importantly, the coefficients B\ and B 2 depend on the orientation of the 
'polarization wave' given by the angle 9. 



ORIENTATION OF THE DOMAIN STRUCTURE 



Consider the domain structure formed close to the paraelectric-ferroelectric transition. Although stability of the 
paraelectric phase is lost with respect to the polarization waves with the value of the k vector given by Eq. (f2T))) and 
arbitrary orientation in the x — y plane, i.e. for any 9, the energy of sinusoidal domain structure depends on 9 and one 
has to find the ones corresponding to the equilibrium domain structure(s). Since we consider here one 'polarization 
wave' only, we study the competition between the stripe-type structures. In Sec lVIII we shall show that the square 
(checkerboard) domain structure is unstable in perovskite crystals that we study here. The checkerboard structure 
could, in principle, be stable or metastable under some conditions on the material constants, but we are not aware of 
any experimental example of this type, so it will be premature to study such a hypothetical case. 

Recall that we discuss the ferroelectric phase transition in a sample with short-circuited electrodes. Then, p = 
in the ferroelectric phase at least not far from the phase transition, and the phase transition into the inhomogeneous 
domain phase occurs at A 2 — 0. The free energy is: 



T/ -ip, s A 2 2 9B 2 (9) 4 
V k(jp,a) = —a A ^rr, — a 



256 



At a fixed 9, the minimum of this free energy is realized for 



W A 2 
9B 2 {9Y 



with the corresponding free energy 



4-' 

~9B 2 (9)' 



(58) 



(59) 



(60) 
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We see that the equilibrium domain structure is realized for the angles 9 which minimize the function B 2 (9) or, 
according to Eq. (|57p . maximize the function / (9). Let us find maxima of this function. It can be written in the form: 



where 



2^ 4Ai/x 

c =m — r~\ — I o.am rr- ( 62 ) 



Ai-A 2 ' (Ai + A 2 + 2n) (Ai - A 2 ) 
One sees that for r > c or 

A 2 +2m>Ai, (63) 
the equilibrium domain structure corresponds to 9 eq = 0, 7r/2 and 

/W roax = /(0) = 1/A 1) 
while in the opposite case 9 eq = it J A, 3tt/4 and 

/(^) max = /W4)= Ai + A 2 2 + 2M - (64) 

Throughout the present paper, we use the data for the material constants of BaTi03 and PbTi03 from Refs&iSr— 
and of Pb(Zr .5Ti .5)O3 from Ref£. We see that for BTO and PTO the condition of Eq. ([63|) is met and, therefore, 
the equilibrium 180° domain structure consists of stripes parallel (perpendicular) to the cubic axes, while the opposite 
inequality applies to PZT and the stripes make 45° with the cubic axes there. For BTO, our conclusion coincides 
with that of Dvorak and Janoveoi^ who defined the equilibrium orientation of the 180° domain walls in BTO far from 
the phase transition. These authors were surprised by their conclusion about a very weak orientational dependence of 
the domain wall energy given that the experimental observations^ showed a clearly preferable orientation, the same 
as suggested by the theory. 

It follows from our results that the weak orientational dependence of the domain structure energy takes place in 
the sinusoidal regime too, and not only for BTO, but for all three perovskites we have made the numerical estimates 
for. Indeed, from Eq. (|60[) one sees that the orientational dependence of the domain structure energy comes from the 
function B 2 (9) . The maximum difference of values Eq. (|5?1) for B2{9) can be used to characterize the anizotropy of 
the domain wall energy: 

AB 2 = B 2max - B 2min = \ q \ 2 |/ (0) - / (tt/4)| = tf ^T^Vf 1 - (65) 

3 3 Ai Ai + A2 + 2/i 

We found AB 2 /B 2 ~ 4 x 10~ 3 for BTO, much smaller anizotropy - 3 x 10~ 4 for PTO, and even smaller one for PZT 
where AB 2 /B 2 ~4x 10 -5 (we have used the parameters listed in the footnote—). Such a weak angular dependence 
of the domain structure energy is in accordance with phase field results of Ref.™ that shows domain walls mainly 
with thermodynamically favorable orientations but also with strong deviations from them. 



VI. LOSS OF STABILITY OF A SINGLE DOMAIN STATE 



It is convenient to study the loss of stability of the single domain state with respect to formation of a domain 
structure using Eq. (|56p. With this, we mean the loss of stability with respect to an arbitrarily small 'polarization 
waves' so that the original single domain state may be, in principle, either stable or metastable. Specifically, in our 
case, when Eq. (|29|) is valid, this state is metastable^. 

A solution of the equations, 

dF/dp = 0, dF/da = 0, (66) 

corresponding to a single domain state (p ^ 0, a = 0) is possible only if A\ < with p\ xtr = —A\/B, where the 
subscript stands for the 'extremum'. This extremum is a minimum (which is relative in our case) if 



d 2 F/d P 2 > 0, d 2 F/da 2 > 0, 



(67) 




FIG. 3: Schematic of the sinusoidal domain structure in (a) BTO and PTO, and (b) PZT. While the stripes are oriented along 
the crystal axes in case (a), in case (b) the stripes are at 45 degrees with respect to cubic axes, the difference being due to an 
opposite sign of elastic anizotropy in those two cases. 
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at the point p = p ex tr, a — given that d 2 F/dpda is evidently zero at this point. The first inequality in Eq. ([67|) 
obviously valid for Ay < 0, while the validity of the second is not immediately evident. 
We find from lEffil): 




= A 2 + 3.Bi (9) p\ xtr =A 2 - 3^i - 3Ai (fli (0) - b) /B, (68) 
From the condition I d 2 F/da 2 ) = 0, we obtain the value of A corresponding to a loss of stability of the single 

V / a=0,p=p cxtr 

domain state with respect to appearance of a polarization wave with a given orientation, A pw (8). It is convenient to 
present it in the form: 

A pw (9) = ^— + Gk 2 + ~ 2GkA (3 (8) , (69) 



£bd + e e l \Ebd + e e l 



where 



f)(0)= ~ g "~ g ' 2/WAl . (70) 

BX 1 +2[q 2 1 -q 2 2 f(9)X 1 } 

The last term in Eg. (1691) is the result of the polarization-strain coupling while the first two present the prior case 
without this coupling^. According to Eq. (j6"T|) . the corresponding single domain state will be (meta) stable at low 
temperatures such that A < min A pw (6) . 

The actual loss of stability of the single domain state corresponds to the minimum of A pw (9) . We have seen in 
Sec|V]that for perovskites the angular dependencies are very weak and we can neglect it putting fX± = 1. Then 

(3 = = g " ~ g ' 2 (71) 

BX 1 +2(qf 1 -qf 2 ) 

Since in BTO, PTO, and PZT q 2 x > q 2 ^, the last term in Eq. (f69"l) is positive and, therefore, the region of 
mefastability of the single domain state in these systems is broader than according to£ and£, in apparent accordance 
within. However, there are serious reservations. First of all, the effect is not very spectacular. Indeed, the factor (3 in 
the last term of Eq.(|69|) is always less than one half, j3 < 1/2, approaching that value when q\ x — q\ 2 — > oo. Therefore, 



And 

A pw < — -, (72) 

Ebd + e e l 



where the r.h.s. corresponds to a very strong strain coupling. Recall that A = —And/ (sbd + e e l), or A\ = 0, 
corresponds to what was calculated in several papers as a 'critical thickness of single- domain ferroelectricity 1 l^ D , 
Eq.pT]). 

To get the opposite limit of a weak strain coupling for A pw (0) , we put q\ Y — q\ 2 — and neglect Gk 2 , as Pertsev 
and Kohlstedlr^ did. We see that 

6ird , , . 4-7TG? , , 

-^ <A - (0)< -^l> (73) 

i.e. because of account for the polarization-strain coupling the value of A pw (0) changes always by less than 1.5 times. 
In usual situation when e&d < e e l this is also the interval of change of the thickness corresponding to absolute loss of 
stability of the single domain state at a fixed temperature. 

The above moderate, less than 50%, range of change is in striking disagreement with a statement by Pertsev and 
Kohlstcdt 10 who claimed a more than an order of magnitude change by nullifying the electrostrictive constants. They 
do not report the details of their procedure, but it is clear from the rest of the paper that their suggestion of putting 
the electrostrictive constants to zero implied changes in the coefficients of the LGD free energy that should have been 
renormalized by the misfit strains. Evidently, it has nothing to do with the effects of the polarization-strain coupling 
omitted in Refs® ,? , since this renormalization is automatically taken into account there, while the effect of misfit strain 
on LGD coefficients was apparently neglected in a gedanken exercise performed in Ref^. 

Specifically, we find that for the perovskites BTO and PTO f3 — 0.4, while in PZT this parameter is 0.1, i.e. four 
times smaller. We see that BTO and PTO are similar and closer to the limit q\ x — q\ 2 — > oo, (3 = 0.5, i.e. the point 
of loss of stability of single domain state is quite close in these materials to the 'critical thickness of single-domain 
ferroelectricity' studied in Refi^ and elsewhere. However, the latter does not have any practical importance because 
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FIG. 4: (color online) Regions of (meta)stability of the single domain and polydomain states in the ferroelectric film as a 
function of temperature T at fixed thickness (top) and as a function of the film thickness I at fixed T (bottom). Upon lowering 
the temperature at a fixed thickness I (top panel), the paraphase gives way to domains that are stable at all temperatures 
T < Td, where Td is below the critical temperature of the bulk ferroelectric transition T c . The single domain (SD) state is 
metastable at low temperatures T < T^f , when a strain coupling is neglected. The strain coupling shifts the boundary of 
metastability towards so-called critical temperature for a single domain state, Tj? D , as marked by the vertical arrows for the 
perovskites in question. The phase behavior of the films as a function of their thickness I at fixed temperature (bottom) is 
qualitatively similar. Very thin films are in a paraelectric phase that is replaced by domains at larger thicknesses I > Id- The 
single domain state is metastable at thicknesses I > ifj^ and becomes suitable for memory applications at even larger (yet to 
be determined) thicknesses when the life time of the metastable state becomes sufficiently large. Strain coupling may extend 
the boundary of the metastability down to the so-called 'critical thickness for SD ferroelectric state' lf D . The films with FE 
memory would be somewhere at larger thicknesses in the metastable region, as marked on the diagrams. 



if one fixes the temperature and reduces the film thickness starting with a monodomain ferroelectric state at low 
temperatures or large thicknesses, that state will give way to domains before the thickness determined by the limit of 
the single domain state stability is reached. The fact of the matter is that the single domain state is metastable, it 
may have a large lifetime at low temperatures and large film thicknesses but this lifetime goes essentially to zero (to 
'atomically' short times) when the above mentioned temperature or thickness are reached. 

Importantly, it follows from Eq. (JHSJ) that if gf x < q 2 2 , the account for the inhomogeneous strains shrinks the 
region of metastability of the single domain state. This shows that contrary to the claim by Pertsev and Kohlstedt 
there is no general physical phenomenon such as stabilization of a single domain state because of inhomogeneous 
strains accompanying formation of domains This may seem surprising, because solids are known to 'dislike' the 
inhomogeneous strains (free energy usually goes up). Moreover, the expectation of Pertsev and Kohlstedt is justified for 
a free-standing him, at least for elastically isotropic solid 13 . But it is not certain for a film on substrate considered both 
by them and in the present work. To explain the physical reason, we recall that the coupling with strain renormalizes 
the coefficients before fourth-order terms in the LGD free energy, in our case we mean the coefficients before p 4 
and p 2 a 2 terms. Then, one has to take into account that the homogeneous strains in the plane of the substrate are 
not possible while inhomogeneous ones are. Both homogeneous and inhomogeneous polarization create homogeneous 
strain but to a different extent see Eqs. (|3"rj]) and (|3"5)l .while inhomogeneous strains are created, of course, by the 



14 



o 

o 



j BTO/SRO/STO \ 




1 2 3 4 5 6 

FE film thickness / (nm) 



FIG. 5: (color online) The phase diagram for BaTiOg/SrRuOs/SrTiOs films in the plane temperature-thickness. The line 
marked Id delineates the para- and domain phases, while the one marked l ms marks the boundary of the metastability regions of 
the single domain state calculated for BTO accounting for strain coupling, Eq. (|69|l . The broken line at larger thicknesses from 
l m s marks the metastable region calculated without accounting for the strain coupling, Eq. (|76[) . The films with FE memory 
would be somewhere at larger thicknesses in the metastable region. 



inhomogeneous polarization only. The out-of-plane and in-plane strains couple with the ferroelectric polarization P z 
by electrostriction terms with different coefficients and the final result of renormalization of the coefficient of p 2 a 2 
term is due to several contributions and it is not clear upfront. It should be obtained by a consistent analysis, as 
it has been done above. No reason is seen to discard the possibility that the inequality q 2 x < q 2 2 can be realized in 
some systems, and one cannot exclude, at least for the moment, the possibility of favoring the multidomain state by 
the polarization-strain coupling. Interestingly enough, this favoring may be very strong: according to Eq. (|69p , the 
increase of the region of absolute instability of single domain state becomes infinite when q\ 2 — q\ x tends to BXi/2 
from below. 

The 'phase diagrams' for the epitaxial FE films on a misfit substrate are plotted in Figs. 0] [5] The boundaries of 
paraelectric phase, domains, and metastable single domain region for BaTiOa/SrRuOs/SrTiOa system are shown in 
the temperature-film thickness (T — I) plane in FigJS] They are found from the conditions that we discussed above 
and write down here for a reference: 



A - 
A -- 
A -- 

where the Landau coefficient A is evaluated 



_ 2Gfc2 _ 4 (^y 2 i (74) 



67Td (>7R\ 

= -TJsd> ( 76 ) 

c e l ms 

for a given temperature T of interest, d/2 = A = 0.8A, e e = 8.45 for 
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SrRu03 electrode, G — 0.3A 2 , and e±_ the dielectric constant [see its definition below Eq. (fl7|) ] in the plane of the 
FE film has been found from the Landau coefficients 5,6 . The last condition corresponds to I = lf^ found without 
accounting for the strain coupling. The critical line T — lf^ in the phase diagram, Fig|5]for BTO with an account for 
strain coupling has been found from Eq. (l5^|) . The arrows on Fig JS] show the evolution of the state at either T = const 
or I = const. 

One should understand that in both illustrations it is implied that the corresponding critical points have physical 
values as solutions to the conditions (??), or, equivalently, Eqs. (f^5j) . (l^71l . (|5Uj) . (|3ip . Consider first the lowering of the 
temperature at a fixed thickness / (Fig. 01 top panel), where the paraphase transforms into domain state below the 
temperature Td that is smaller than the critical temperature of the bulk ferroelectric transition T c . We see that the 
single domain (SD) state would be metastable at low temperatures T < T^f in the region overlapping with the 
domain state. Note that the Tjjj^ plotted in Figj5]is found without accounting for the strain coupling. The strain 
coupling then shifts the boundary of metastability towards the so-called critical temperature for a single domain state 
Tjf D thus broadening the range of metastability of the SD state, as shown in Fig|4j The phase behavior of the films 
as a function of their thickness I at fixed temperature (Fig. 01 bottom panel) is qualitatively similar. Very thin films 
remain in a paraelectric phase that is replaced by the domains at larger thicknesses I > Id- We see that both T C 5D 
and if D are actually unreachable in the present case, since the system may get to those points only by moving from 
the paraphase down (right to left on the phase diagram, FigJS]), but such transitions are preempted by the domain 
instability that sets in first. The single domain state is metastable at thicknesses I > lf^, and becomes suitable for 
memory applications at even larger (yet to be determined) thicknesses where its life time becomes sufficiently long. 



VII. INSTABILITY OF THE CHECKERBOARD DOMAIN STRUCTURE 



In the previous Sections, we have assumed that the domain structure is stripe-like by taking into account only 
one 'polarization wave'. This and other possibilities have been studied by Chensky and Tarasenko^ who considered 
the uniaxial ferroelectric isotropic in the x — y plane. Along with the stripe-like structure they discussed also the 
checkerboard and the hexagonal domain structures. The latter can be realized in the presence of an external field only 
which is not discussed in this paper. However, a checkerboard structure should be analyzed as an alternative to the 
stripe structure. In Ref^, the authors stated that the checkerboard structure never realizes, although, surprisingly, 
there is no proof of this statement. In this Section, we shall show that this structure is indeed unstable for the 
isotropic case treated in Ref£ and then show that this conclusion holds also when one explicitly takes into account 
the polarization- strain interaction, apart from mere renormalization of the LGD coefficients by the misfit strains. 

Once again, we consider a short-circuited sample, i.e. the ferroelectric polarization is described by: 

P z = a\ cos k\r cos qz + a 2 cos k 2 r cos qz, (77) 

where k\ and k 2 are two noncollinear vectors whose modulus is given by Eq. (|20p and whose (mutually orthogonal) 
directions remain unspecified for a moment, Fig. [5J 



A. Checkerboard domains without elastic coupling (gn = q\2 = 0) 



In this case the solution for the fields is [cf. Eq. (fT0)) ] 

E z = E^ 1 cos kir cos qz + E z 2 cos k 2 r cos qz , (78) 
E X /y = E^]y sinfcir singz + 23* 2 sin k 2 r sinqz, (79) 

with Eg. (1181) still applicable to spatial harmonics (as follows from linearity of Maxwell equations) and the equation 
of state for the fundamental harmonics is the same as Eq.([Tj 



A 2 a 1 +[BP*] cc = E*\ (80) 

where one retains the terms [.BP.^] oc cosfcircosgz (symmetry dictates the analogous expressions for 02). In the 
above equation, 

-P? = (ai cos kir + a 2 cos k 2 r) 3 cos 3 qz — — (a? + 2aia 2 ) cos k\r cos qz + — (a 2 + 2a 2 af) cos k 2 r cos qz + . . . , (81) 

16 16 

so that we obtain for the fundamental harmonic the following equations of state: 

9B 

A 2 a x + — (a\ + 2 ai al) = (82) 
lb 
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FIG. 6: Schematic of the checkerboard domain structure. It is absolutely unstable in case of BaTiOs, PbTi03, and 
Pb(Zro.5Tio.5)03 typical perovskite ferroelectrics. 



and the analogous equation for a 2 . Since these equations are obtained from extremum of the free energy, dF/dai^) = 
0, we again restore the full free energy, accounting for the symmetric contribution by the a 2 harmonic: 

V-'F = ^(af + a%) + ^B {a\ + af) + ^Balal (83) 

The equations of state have the checkerboard solution: 

a\ = a\ = -16A 2 /27B. (84) 
Checking what type of extremum for the free energy is this solution, 

d 2 F d 2 F f d 2 F \ 2 _ l_ 2 
da\ X da\ \daxdai) ~ 12 2 ' 

we see that the checkerboard solution is the maximum of the free energy for some directions in the 01,02 plane and 
is absolutely unstable. 



B. Checkerboard domains with elastic coupling 



We have seen above that the elastic coupling renormalizes the fourth order coefficients in formulas like Eq. (|83|) 
reducing them by some amounts which are different for different coefficients. Thus, instead of Eq. (|53"|) . we will have: 

V-'F = ^(aj +4) + ^ (B 21 af + B 22 ai) + ^B 3 afal (85) 

where B 2 \ and B 22 are given by Eq. ([57[) for the corresponding angles and B3 is a new coefficient which depends on the 
both angles and which can be, in principle, both positive and negative. Both from Eq. J57[) and the cubic symmetry 
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one realizes that B 2 \ = B 22 = B 2 - For what follows, it is important to mention that when B 3 is negative it cannot be 
of large absolute value, otherwise there will be directions in the (a 1; 02) plane along which the free energy diminishes 
without limits at large values of a\,a 2 what means a global instability of the system. By putting a\ = 02, one sees 
from Eq. (|55|) that to avoid this instability the condition 



B 2 + 2B 3 > 



(86) 



should be fulfilled. Another evident condition of the global stability is B2 > 0. 
The checkerboard solution is: 



-I6A2 



' 9 (B 2 + 2B 3 ) ' 

To analyze stability of this solution, we calculate the second derivatives 



(87) 



d 2 F d 2 F 



daf 



da 2 



A 
4 



27 



Bo 



-A 

2 A B 2 + 2B 3 7 



d 2 F 9_ B 
da\da 2 16 



±- 



A 2 B 3 
Bo + 2Bo 



we then find the discriminant 



da 2 )\da 2 ) {darfazj A 2 B 2 +2B 3 1 ' 

Our further study is aimed at finding out if and when the condition of positiveness of Z, i.e.B 2 > 2B 3 , is compatible 
with the two conditions of the global stability mentioned above. Thus we need formulas for the coefficients B 2 and 
/>'•,. 

Turning to taking into account explicitly the polarization-strain coupling, we recall that in our approximation of 
sufficiently thick film the only source of the elastic strains is P 2 . This function contains now a cross term stemming 
from 



P 2 — [a\ cos kir cos qz + a 2 cos k 2 r cos qz) (89) 

1 ~\~ cos *2,c[z 

= [a 2 + a 2 , + 2a\a 2 (cosp + r + cosp~r) + a\ cos 2k\r + a 2 cos 2fc2r] , (90) 

where = k\ ± k 2 . Naturally, the components of the strain tensor will have terms depending on cosp + r, cosp _ r:. 
We will have for u xx 

u X x — "ix + U S cos + u xt cosp + r + cosp~r + u xx cos p + r cos 2qz 

+u%~ cosp^r cos 2qz + uf} xx cos 2fcir + u^l x cos 2fc 2 r + f^i^L cos 2k ± r + u^ xx cos 2k 2 rj cos 2qz, 

and the analogous equations for u yy and u zz . 

From our previous experience, it becomes immediately clear that u ^(yy) = u xx(yy)> smce the equations for those 
components do not depend on z and we can now drop the indices p and q from the corresponding terms, leaving only 
u + 'T v Also, due to the same reason as above, u^°\ <, = , = 0. Then, one can write: 

xx(yy) ' ' xx(yy) xx{yy) ' 

u X x = (uxx cosp + r + u xx cosp~r) (1 + cos2gz) 

+Uxl x cos 2k\r + u^xx cos 2k 2 r + (u^lx cos 2k\T + u^xx cos 2fc2»") cos 2qz. 

and a similar equation for u yy . The 'diagonal' terms for the first (second) feu2) harmonic are 



(3) (3) _ (4) (4) _ _ a l(2)f 

U l(2),xx + U l(2), yy ~ U l(2),xx + U 1(2), TO ~ 912 . 
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All the cross terms satisfy 



9 2 u± d 2 u± d 2 u ± dP 2± 

^^ + (^+^^+^^ + 912^=0, (91) 
ox z oyox oy z ox 



d 2 u± „ 8 2 u± d 2 u± dP 2± 

X ^ + ^ + ^dytx + ^ + ^^y~ = °' ™ 

where P 2± — \a\a 2 cosp ± r, in components u^,^ oc sinp ± r, d 2 uf/dx 2 — — (pf) 2 uf etc. 

^wt 2 ut + (Aa + n)ptPyUy + fiPy 2 u± + q 12 p~t\a\ai = 0, 



(93) 



^lpt 2 Uy + (A 2 + Ll)ptPyut + [-ipt 2u y + qi2ptn a l a 2 = 0- (94) 



J y "-x 1 r^i "•y 1 Hi-^fy ^ 

Then, 



^xx + u yy Pa; w x + Py u y 



-?i2^/(^). (95) 

.de, as above, that only -u^and u£s are non-zerc 

phase, we obtain 



For u zz we conclude, as above, that only u z °} and u£~) are non-zero and following the same reasoning as for the stripe 



u ( z °J=uW = -^-(a 2 + a 2 ). 

Having solved the elastic problem, we are now in a position to write down the constituent equations containing a\ 
and a 2 only. To this end, we write two equations for a\ and a 2 analogous to Eg. (1211) but this time [. . .] would mean 
the proportionality to cos k\r cos qz or cos k 2 r cos qz, respectively. Since both equations have the same structure, we 
will discuss that for a\ only and, for the sake of brevity, we will mention only the terms containing a 2 , the other terms 
are the same as for the one-sinusoid case discussed above. 

For clarity sake, we repeat Eq. ([2l"T) with a minor change for the present case: 

A 2 a x + [BP3 + 2q n P z u zz + 2q 12 P z (u xx + u vv )} cc = 0. (96) 

It is straightforward to find that the a2-containing term stemming from [P z ] is 9aia 2 /8. From Eg. (1^51) . one sees 
that [P z u zz ] = a± (uzz + \ u ^tz ) , recall that now we consider the case p — 0, and the contribution of this term is 



3gn 2 
~8A7 aia2 - 

Now 

3 



[P z {uxx + u yv )] cc = y (u[ 3) xx + uf^ + y (u[% + ufl^j + -a 2 (u+ x + u+ y + u xx + u yy ) 



and contribution of this term to the equation of state is 



-* 12 ^[/(0 + )+/(n] 



Finally, the constituent equation for a\ takes the form: 



where we have introduced 



A 2ai + ^a\B 2 {Ox) + ^a ia 2 2 B 3 (9+ ,0~) = 0, 



B 3 (0\6-) =B-\ q f i - \q\ 2 [f(9 + ) + f(9-)] (97) 
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Similarly to the case of one sinusoid, we recover the free energy: 

F = ^ K + al) + ^B 2 (9,) 4 + lfi 2 (9 2 ) 4 + ^B 3 {6+,9~) a\al (98) 



Recall that the square symmetry suggests that B 2 (9\) = B 2 (9 2 ). and B 2 (9) is given by Eq. (|57|) 

Turning to examining the sign of the discriminant Z, we should mention that according to Eqs. ([57| and (|9"T|) : 



B 2 - 2B 3 = —B + ^q\ 2 



f{9 + ) + f{9-)- l -f{9 1 ) 



One sees that the maxima of Z correspond to maxima of f{9 + ) and f(9 ) [note that f(9^ lax ) — /(# max ) because of 
the cubic symmetry], which are, automatically, the minima of / (#1) as we have seen in SecJV] Then, 



[ft- = --B+~«? 2 



Using the values of f(9 max ) and / (0 m i n ) found in Sec|V] we find that if X 2 + 2\i > Ai, 

4 , / 2 1 \ 4 , 2 (A 2 + 2u) + Ai 



4 a / 4 1\ 4 2 7Ai - A 2 - 2/i 



and if A 2 + 2\i < Ai, 

- 2B 3]max = —b + ^ Ai+A ; +2 ^ - = -5 + 3& 2A ;- (Ai+ - 2+ -^ doo) 

To prove that the checkerboard structure can be stable, at least in principle, with respect to small fluctua- 
tions, we should demonstrate that the positiveness of [B 2 — 2£?3] max is compatible with the conditions B 2 > and 
[B 2 — 2i?3] min > which guarantee the global stability of the system. We do not intend to perform an exhaustive 
analysis but want only to demonstrate that this is possible under certain conditions, unlike in the case without the 
elastic coupling. As an example, we consider a system with a weak elastic anisotropy which is valid for the perovskitcs, 
i.e. we shall assume X 2 + 2fj, ~ Ai, and q%\ = 0. Both Eq. ([M|) and Eq. (|100l) then give 

[ft -2B a ] mKX ~-B + ?&, (101) 

and give for the positiveness of [B 2 — 2i?3] max the same condition:^ > BXx/2 while the condition B 2 > now reads 
<7i 2 < 3SAi/2. One sees that for a nearly elastically isotropic ferroelectric with qu = 0, the checkerboard structure is 
at least metastable if 

3SAi/2 > q\ 2 > BXi/2 

Of course, the above set of the material coefficients looks fairly exotic but it is just an example aimed at nothing more 
but demonstration that the checkerboard domain structures are permitted due to the elastic coupling when certain 
conditions on the material coefficients are met. 

In case of real perovskite films the checkerboard structure is not stable. To see this, we can rewrite Eq. (|101[) in 
the form: 

[ft - 2B 3 ] max * -S - ^ _ ^ < o 

Indeed, we have already mentioned above that for the perovskites q\ x > q\ 2 . Also, B > there. Therefore, in the 
perovskites the checkerboard domain structure is absolutely unstable 



VIII. CONCLUSIONS 



With the use of the Landau-Ginzburg-Devonshirc theory, we have studied the effects of polarization-strain coupling 
when defining the character of equilibrium domain structures and the limits of absolute instability of a single domain 
state in thin films of cubic ferroelectric films on a misfit substrate. On the compressive substrate, the cubic ferroelectric 
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behaves substantially as a uniaxial ferroelectric with the polar axis perpendicular to the film. The film is sandwiched 
between the electrodes that do not provide a perfect screening of the depolarizing field because of the finite Thomas- 
Fermi screening. length. Such a system is exemplified by (100) BaTiOs/SrRuOs/SrTiOs film and similar perovskite 
structures. Quantitative results have been obtained for BaTiC>3, PbTiC>3, and Pb(Zro. 5^0.5)03. We have found 
that close to the paraelectric-ferroelectric phase transition or at the film thicknesses close to the minimal thickness 
compatible with the ferroelectricity, the equilibrium domain structure in perovskites is the stripe-wise one with 
the stripes parallel (perpendicular) to the cubic axes in BaTi03, PbTiC>3, while running at 45° to cubic axes in 
Pb(Zro. 5^0.5)03. The energy of the domain structure depends very weakly on the stripe orientation, the maximum 
change proves to be well below 1% in all three cases. We found that because of the polarization-strain coupling 
a competing checkerboard domain structure may, at least in principle, be equilibrium or metastable when certain 
conditions on the material constants are met, but we are not aware of any material system with such conditions. 
The limit of absolute instability of the single domain state changes due to the polarization-strain coupling. Thus, the 
interval where the absolute instability is absent, meaning a metastability in the cases at hand, widens in perovskites 
in agreement with the earlier conclusion by Pertsev and Kohlstedt 10 . However, this effect is much smaller than 
that claimed by them. Increase of the metastability range is substantial in BaTiC>3 and PbTiC>3 where the absolute 
instability limit becomes close to what is often called the "critical thickness for ferroelectricity" Zf D ' £ FigJSJ but 
without accounting for the domain formation, which sets in first and prevents the system from ever reaching this 
point. The effect is much smaller in Pb(Zro. 5^0.5)03. We have found also that the polarization-strain coupling can 
lead to narrowing of the region of relative stability of the single domain state under certain conditions on the material 
constants, but we are not aware of an experimental. realization of these conditions. 

APL has been partially supported by Ministry of Science and Education of Russian Federation (State Contract # 
02.740.11.5156). 
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